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Abstract 

In the previous paper (Part 1), we have verified that the SK assumption on the direction does not hold in 
the analysis of neutrino events occurred inside the SK detector, which is the cornerstone for their analysis 
of zenith angle distributions of neutrino events. Based on the correlation between L v and (Figures 16 
to 18 in Part 1) and the correlation between E v and (Figure 19 in Part 1), we have made four possible 
L/E analyses, namely L v /E Vl L u /E^ : L^/E^ and L^/E^. Among four kinds of L/E analyses, we have 
shown that only L v /E v analysis can give the signature of maximum oscillations clearly, not only the first 
maximum oscillation but also the second and third maximum oscillation and etc., as they should be, while 
the L^/E^ analysis which are really done by Super-Kamiokande Collaboration cannot give any maximum 
oscillation at all. It is thus concluded from those results that the experiments with the use of the cosmic-ray 
beam for neutrino oscillation, such as Super-Kamiokande type experiment, are unable to lead the maximum 
oscillation from their L/E analysis, because the incident neutrino cannot be observed due to its neutrality. 
Therefore, we would suggest Super-Kamiokande Collaboration to re-analyze the zenith angle distribution of 
the neutrino events which occur inside the detector carefully, since L v and L^ are alternative expressions of 
the cosine of the zenith angle for the incident neutrino and that for the emitted muon, respectively. 
PACS: 13.15.+g, 14.60.-z 
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1. Introduction 

The specification of the oscillation parameters 
for neutrino oscillation is entirely based on the sur- 
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vival probability for a given flavor (Eq.(l)) in which 
two physical quantities to be measured, namely, the 
directions of incident neutrinos and their energies 
are included. However, these two quantities cannot 
be measured directly due to their neutrality and 



ures with colors are strongly impressive compare with those 
with monochrome. In the figures with colors, we classify 
neutrino events by blue points and aniti-neutrino events by 
orange ones. 
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Super-Kamiokande Collaboration are forced to in- 
troduce the SK assumption on the direction. 

As shown in Figures 11 to 13 and/or Figures 16 
to 18 of the preceding paper [![, we have shown 
that the SK assumption on the direction requiring 
that the directions of the incident neutrinos are the 
same as those of the emitted muons docs not hold 
in the case of the neutrino events with the high- 
est quality, namely the single ring muon events due 
to quasi-elastic scattering(QEL) among Fully Con- 
tained Events. 

Also, in Figure 19 of the same paper, we have 
shown that the energies of the incident neutrinos 
cannot be determined uniquely from those of the 
emitted muons. 

Compared Figures 11 to 18 with Figures 19 in the 
preceding paper l], it is easily understood that the 
approximated E v by Eq.(9) in [l[ does not bring a 
fatal error into the survival probability for a given 
flavor (Eq.(l) in the present paper) in spite of the 
unsuitable theoretical treatment, while the SK as- 
sumption on the direction introduces the fatal error 
into the L/E analysis. 

In the present paper, we examine how does the 
variable L/E in the survival probability for a given 
flavor influences over the results around L/E anal- 
yses. In our computer numerical experiments, it 
is possible to analyze four different kinds of L/E 
distribution, namely, _L„/_E„, L^/E^, L^/E^ and 
L^/E v distributions, where L„, L M , E v and de- 
note the flight length of neutrino, the corresponding 
flight length of emitted muon, the incident neutrino 
energy and the emitted muon energy in QEL, re- 
spectively. 
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Figure 1: L v /E„ distribution without oscillation for 1489.2 
live days (one SK live day). 
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Fi gure 2: Ly/Ei, distribution without oscillation for 37230 
live days (25 SK live days). 



2. L/E Distributions in Our Computer Nu- 
merical Experiment 

Here, we explain the procedure of our computer 
numerical experiment roughly (See details in the 
Appendices in the preceding paper [l[). At first we 
construct a hypothetical SK detector in the com- 
puter. We randomly sample a neutrino event with 
both certain energy and zenith angle from neutrino 
spectra at the opposite side of the Earth, imaging 
the injection of the neutrino concerned into the de- 
tector. We pursue the neutrino concerned up to 
the inside detector where the neutrino interaction 
is expected^- In the interaction (QEL) occurred in- 



2 Exactly speaking, instead of the description in the text, 
we adopt to sample the neutrino events from the neutrino 
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Figure 3: Survival probability of P(uu — > f M ) as a function 
of L v / E u under the neutrino oscillation parameters obtained 
by Super-Kamiokande Collaboration. A, B and C represent 
the first, the second and the third maximum oscillation, re- 
spectively . 
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side the detector, we "measure" the muon energy 
from the random sampling of Q 2 for the neutrino 
concerned and we pursue the muon concerned up 
to the end of the detector, taking into account of 
all physical processes due to the muon and judging 
whether the muon concerned belongs to Fully Con- 
tained Event or Partially Contained Events. There, 
we adopt Fully Contained Events only. As the result 
of a series of these procedure, we are able to know 
a series of a pair of the neutrino with the known 
primary energy, E v , and the zenith angle, cos9 v or 
L„, and the produced muon with the energy E^ 
and the zenith angle, cosO^ or L^. In our computer 
numerical experiments, every physical process is 
treated stochastically and physical results are thus 
obtained, taking account of the stochastic charac- 
ters inherent in their processes exactly. Namely in 
this sense, there is one-to-one correspondence be- 
tween "measured" neutrinos and "measured" their 
daughters muons, while one cannot generally spec- 
ify the parent neutrinos from the measured muons 
in the real experiment. 

Our computer numerical experiments are carried 
out in the unit of 1489.2 days. The live days of 
1489.2 is the total live days for the analysis of the 
neutrino events generated inside the detector used 
by Super-Kamiokande Collaboration 2] . Hereafter, 
we call 1489.2 live days as one SK live day. We 
repeat one SK live day experiment as much as 25 
times, namely, the total live days for our computer 
numerical experiments is 37230 live days (25 SK 
live days). 

2.1. L v /E v distribution 

2.1.1. For null oscillation 

In Figure [TJ we show L v jE v distribution with- 
out oscillation for one experiment (1489.2 live days) 
among twenty five computer numerical experiments 
(25 SK live days). In those numerical experiments, 
there are statistical uncertainties only which are 
due to both the stochastic character in the phys- 
ical processes concerned and the geometry of the 
detector. Therefore we add the standard devia- 
tion as for the statistical uncertainty around their 
average in the forthcoming graphs, if neccessary. 
In Figure [3J we show the statistical uncertainty, 
the standard deviations around their average val- 
ues through twenty five experiments. Similarly for 



interaction spectra inside the detector, which is mathemati- 
cally equivalent to the procedure described in the text. 



other possible combinations of L and E ( L^/E^, 
L^/E^ and L^/Eu) for 37230 live days (25 SK live 
days) we did so. 

In Figures Q] and [2j both distributions show 
the sinusoidal-like character for L v jE v distribution, 
namely, the appearance of the top and the bottom, 
even for null oscillation. In this case, it should 
be noticed that their distributions are expressed 
in a logarithmic scale. The uneven histograms in 
Figure [TJ comparing with those in Figure [2j show 
that the statistics of Figure [TJ is not enough com- 
pared with that of Figure [2] Roughly speaking, 
smaller parts of L v /E v correspond to the contri- 
bution from downward neutrinos, larger parts of 
L v / E v correspond to those from upward neutrinos 
and LyfEy near the minimum correspond to the 
horizontal neutrinos, although the real situation is 
more complicated, because the backscattering ef- 
fect in QEL as well as the azimuthal angle effect in 
QEL could not be neglected as shown in the pre- 
ceding paper[l|. From Figured! we understand that 
the bottom around 70 km/GeV corresponds to the 
contribution mainly from the horizontal-like direc- 
tion and has no relation with neutrino oscillation in 
any sense. 

2.1.2. For oscillation (SK oscillation parameters) 
The survival probability of a given flavor, such as 
v,j,, is given by 

P(v M -> Vn) = l-sin 2 29-sin 2 (1.27Am 2 L„/E l/ ). (1) 

Then, for maximum oscillations under SK neutrino 
oscillation parameters [2|, we have 

1.27 Am 2 I„/£v = (2n + 1) x |, (2) 

where sin 2 29 = 1.0 and Am 2 = 2.4 x 10 _3 eV 2 . 
From Eq.(l), we have the following values of L v jE v 
for maximum oscillations. 

L v /E v = 515km/GeV forn = Q (3 - 1) 
= 1540km/GeV for n=l (3 - 2) 
= 2575km/GeV for n = 2 (3 - 3) 
and so on. 

In Figure [3j we give the survival probability 
P( u fj, ~ * ^v) as a function of L v /E v under the 
neutrino oscillation parameters obtained by Super- 
Kamiokande Collaboration. In cosmic ray experi- 
ments, the energy spectrum of the incident neutri- 
nos is convoluted into the survival probability. 
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Fi gure 4: L^/Ei/ distribution with, oscillation for 1489.2 live 
days (one SK live day), sample No.l. 
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Fi gure 5: EujEu distribution with oscillation for 1489.2 live 
days (one SK live day), sample No. 2. 
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Figure 6: L v /E v distribution with oscillation for 14892 live 
days (10 SK live days). 
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Figure 7: Ly/E^ distribution with and without oscillation 
for 14892 live days (10 SK live days). 
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Figure 8: L v /E v distribution with standard deviations with 
oscillation for 37230 live days (25 SK live days). 



In Figure HI we give one example of our L u /E u 
distribution for one SK live day (1489.2 live days)0 
among twenty five sets of the computer numerical 
experiments in the unit of one SK live day. In Fig- 
ure 03 we give another example for one SK live 
day. Arrows A, B and C represent locations for 
the first, the second and the third maximum oscil- 
lation which are given in Eqs. (3-1), (3-2) and (3- 
3), respectively. By the definition of our computer 
numerical experiments, there are no experimental 
error bars in L v jE v distributions in Figures 0] and 
® 

In Figure [SJ we show the L v jE v distribution for 
14892 live days (10 SK live days). Compared Fig- 
ure E] with Figures H] and it is clear that L v jE v 
distribution in Figure [6] becomes smoother due to 
larger statistics. 

In Figure[jj we give L v jE v distribution for 14892 
live days (10 SK live days) in a linear scale to- 
gether with the corresponding one without oscil- 
lation. The L/E distribution with oscillation in 
Figure [7] is a representation in linear scale and it is 
the same as that in Figure |6] in logarithmic scale. 

However, it is clear in Figure [7] that magnitudes 
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of freqency indicated by arrows A, B and C are 
almost sam^l- It should be noticed from Fgure [7] 
that L v jE v distribution without oscillation repre- 
sents the envelop of the corresponding distribution 
with oscillation. Namely, L v / E v distribution with 
oscillation is equivalent to the L v /E v distribution 
without oscillation multiplied by the survival prob- 
ability for a given flavor (Eq.(6)). 

We have repeated the computer numerical exper- 
iment for one SK live day as much as twenty five 
times independently, in both cases with oscillation 
and without oscillation. In Figure [FJ we can add 
the statistical uncertainty (standard deviation in 
this case) around their average, because every one 
SK live day experiment among twenty five sets of 
the experiments fluctuates one by one due to their 
stochastic character in their physical processes and 
geometrical conditions of the detectors concerned. 

In order to make the image of the maximum oscil- 
lations in L v jE v distributions clearer, we show the 
correlations between L v and E v in Figures 1^1 and [TU1 
which correspond to Figures 0] and [Bl respectively. 
In Figure [5] for one SK live day, we can observe va- 
cant regions for the events concerned assigned by 
A, B and C. In Figure [10] for ten SK live days, the 
existence of the vacant regions for the events con- 
cerned becomes clearer due to larger statistics. 

We give ratios of {L v /E v ) distribution with os- 
cillation to that without oscillation for 1489.2 live 
days (one SK live day) in Figure Qj] and for 14892 
live days (10 SK live days) in Figure [121 respec- 
tively. 

The situation shown in Figures 2] to [TU] shows 
definitely that our computer numerical experiments 
are carried out exactly from the view point of the 
stochastic treatment to the matter, if neutrino oscil- 
lation parameters obtained by Super-Kamiokande 
are correct. 

2.2. L v /E^ distribution 
2.2.1. For null oscillation 

In Figure [T51 we give L u /E^ distribution with- 
out oscillation for 37230 live days (25 SK live days) 
of Super-Kamiokande Experiment to consider the 
statistical fluctuation effect as precisely as possible. 



s Supcrkamiokande collaboration mention that the second 
or more higher maximum oscillations are not observed due 
to their experimental conditional. However, according to 
our results, they could have observed the second well and 
even higher maximum oscillations, if they could have really 
occured like the first maximum oscillation. See arrows A, B 
and C in Fgure [7] 
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Figure 9: Correlation diagram between L„ and E v with os- 
cillation for 1489.2 live days (one SK live day). 
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Figure 10: Correlation diagram between L„ and E v with 
oscillation for 14892 live days (10 SK live days). 
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Figure 11: Ratios of (L v / E u ) osc / {L v / E v ) null for 1489.2 live 
days (one SK live day). 
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Figure 12: Ratios of (L v / E v ) osc /{L v / E v ) null for 14892 live 
days (10 SK live days). 
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Figure 13: Ly/E^ distribution without oscillation for 37230 
days (25 SK live days). 
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Figure 14: Lu/E^ distribution with oscillation for 37230 
days (25 SK live days). 
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Figure 15: L v /E^ distribution with and without oscillation 
for 14892 days (10 SK live days). 



It is seen from the comparison of Figure [T3] with 
Figure [5] for L u /E^ distribution that there is no 
appreciable difference between them and it denotes 
that the transform from to E v (see Figures 19 
in [ljl) does not cause any appreciable change in 
L/E distribution. In other words, the appreciable 
change are caused by the transform of to L v (see 
Figures 11 to 13 ,and/or Figures 16 to 18 in [![). 

2.2.2. For oscillation (SK oscillation parameters) 

Here, we compare Figure IT4l L u /E^ distribution 
for 37230 live days (25 SK live days), with Fig- 
ure [51 corresponding L v jE v distribution. Com- 
bined Figure [T3] with Figure [HI we give L„/£^ dis- 
tributions with and without oscillation are given in 
a linear scale in Figure 1151 Being different from 
Figure [7] for L v /E v distribution, L^/E^ distribu- 
tion with oscillation in Figure [15] begins to make 
the maximum oscillation pattern collapse after the 
first maximum oscillation. This fact corresponds 
to the situation that in L v / E^ distributions, the 
transform of E^ from E v makes it difficult to form 
the " envelope- like" relation between L v jE^ distri- 
butions with and without oscillation after the first 
maximum oscillation. It is clear from the compar- 
ison of these figures that L v /E^ distribution coin- 
cides almost with L v /E v distribution around the 
first maximum oscillation, but the former become 
to deviate from the latter after the second maxi- 
mum oscillation. It reflects the correlation between 
E^ and E v (see Figures 19 in [lj). The situa- 
tion that the first maximum oscillation can be " ob- 
served" is understandable from the existence of the 
vacant region indicated by arrow A in Figure 1161 
too. However, it is needless to say that both L v jE v 
and LyjE^ distributions cannot be observed be- 
cause of the neutrality of L v . 
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Figure 16: Correlation diagram between L v and E^ with 
oscillation for 14892 days (10 SK live days). 



2.3. L^/E^ distribution 

The physical quantities measured by Super- 
Kamiokande Collaboration are and E^ , but nei- 
ther L v and E v . In this sense, we carefully examine 
the validity of the survival probability for a given 
flavor whose variables are and E^ , but neither 
L v and E v . In other words, we examine whether we 
can find the maximum oscillations on L^/ E^ dis- 
tribution or not, because the existence of the max- 
imum oscillation in L^/ E^ distribution is exactly 
the same as the existence of the survival probabil- 
ity for a given flavor whose variable are L M and E^ 
and vice-versa. 
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Figure 17: L M /_E M distribution without oscillation for 1489.2 
live days (one SK live day). 
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Figure 18: L^/E^ distribution without oscillation for 37230 
live days (25 SK live days). 
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Figure 19: L^/E^ distribution with oscillation for 1489.2 
live days (one SK live day). 
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Figure 20: L^/E^ distribution with oscillation for 37230 live 
days (25 SK live days). 



2.3.1. For null oscillation 

In Figure [T71 we give one sample for one SK live 
day (1489.2 live days) from the totally 37230 live 
days (25 SK live days) events, each of which has 
1489.2 live days. FigurelTSlshows the average distri- 
bution accompanied by the statistical uncertainty 
bar (not experimental error bar). It is clear from 
these figures that the existence of the dip or bot- 
tom, namely the sinusoidal character, means the 
contribution merely from horizontal-like contribu- 
tion, having no relation with any neutrino oscilla- 
tion character, as they must be. 

2.3.2. For oscillation (SK oscillation parameters) 

In Figures UHl and [201 we give the L^/E^ distri- 
butions with oscillation for 1489.2 live days (one SK 
live day) and 37230 live days (25 SK live days), re- 
spectively. In Figure fT9l we may observe the uneven 
histogram, something like curious bottoms coming 
from neutrino oscillation. However, in Figure [201 
where the statistics is 25 times as much as that 
of Figure 1191 the histogram becomes smoother and 
such bottoms disappear, which turns out finally for 
the bottoms to be pseudo. 

In order to examine the existence of the maxi- 
mum oscillations in L/E distribution, it is better 
to express L/E distribution in a linear scale, but 
not in a logarithmic scale, as shown in Figure [7] for 
Lu/Ey distribution. In Figure we give L^/E^ 
distribution with and without oscillation for 10 SK 
live days (14892 live days). It is clear from the com- 
parison of Figure [5T] with Figure [7] for L v /E v dis- 
tribution that we cannot find any maximum oscil- 
lation in L^/E^ distribution with oscillation. Also, 
we find the L v /E v distribution without oscillation 
forms an envelop of the corresponding one with os- 
cillation in Figure [7J as it must be, while we cannot 
find such a relation on L^/E^ distribution in Fig- 
ure [2U This denotes that L^/E^ cannot be the 
variables of the survival probability for a given fla- 
vor. In order to confirm the lack of the maximum 
oscillations in L^/E^ distribution, we give a cor- 
relation diagram between L^ and E^ for one SK 
live days in Figure and that for 10 SK live days 
in Figure 1231 respectively. If the maximum oscil- 
lations really exist in L^/E^ distribution, then we 
can expect to find the vacant regions for L M and 
E^ diagrams in Figures |2"21 and |2"31 as shown clearly 
in L v and ^diagrams of Figures [21 and [TOJ How- 
ever, we cannot find anything like vacant regions in 
Figures [22] and [23] at all. 
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Figure 21: L^/E^ distribution with and without oscillation 
for 14892 live days (10 SK live days). 
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Figure 22: Correlation diagram between and with 
oscillation for 1489.2 live days (one SK live day). 
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Figure 23: Correlation diagram between L M and E^ with 
oscillation for 14892 live days (10 SK live days). 



In Figure 1241 we show one example of 
(L^/E^osc/ {L^/E^nuu for one SK live day among 
all possible sets of ratios. We may find pseudo dips 
in the figure. In Figure[25l we give those for ten SK 
live days, whose statistics is larger than that of Fig- 
ure[M]by 10 times, such pseudo dips disappear here. 
Thus the histogram becomes a rather decreasing 
function of L^/E^ in Figure [2SJ When we further 
make statistics higher, the survival probability for 
L^/E^ distribution should be a monotonously de- 
creasing function of L^/ E tll whithout showing any 
characteristics of the maximum oscillation, which 
is in remarkable contrast to Figures Qj] or [T^] for 
L v / E v distribution. 

Summarized Figures[in]to[25l we say that L^/E^ 
distribution cannot give the maximum oscillations 
in any sense. This denotes that L^/E^ distribu- 
tions are not constructed based on the survival 
probability for a given flavor which is the funda- 
mental principle for neutrino oscillation. 



2.1 



L^jEy distribution 



Instead of analyzing L^/E^ distribution, Super- 
Kamiokande Collaboration have analyzed L^/E^ 
distribution where E v is approximated as the poly- 
nomial of E v (See,Eq.(7) in the preceding paper[l|). 
Consequently, we examine the L^/ E v distribution. 



2.4-1- For null oscillation 

In Figuresl2"B1andl2"Tl we give L^/E^ distributions 
without oscillation for 1489.2 live days (one SK live 
day) and 37230 live days (25 SK live days), respec- 
tively. Comparing Figure [5^] with Figure [571 the 
larger statistics makes the distribution smoother. 
Also, there is a sinusoidal-like bottom expressed in 
a logarithmic scale which has no relation with neu- 
trino oscillation. 



2-4-2. For oscillation (SK oscillation parameters) 
In Figures |2"51 and |2"91 we give the L^/E^ distri- 
bution with oscillation for 1489.2 live days (one SK 
live day) and 37230 live days (25 SK live days), re- 
spectively. In Figure [28l we may find something like 
a bottom near ^200 (km/GeV). However, such the 
dip disappears, by making the statistics larger as 
shown in Figure [2"§l At any rate, we cannot find any 
indication on the maximum oscillation from these 
figures. 
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Figure 24: The ratio of (L M /S fl ) osc /(L M /B M ) ntlii for 1489.2 
live days (one SK live day). 



Figure 26: L M /i?„ distribution without oscillation for 1489.2 
live days (one SK live day). 



L/E distribution for FC events 

14892 days Oscillation Am 2 = 2. 4x 10" 3 eV : 







Lu/Eft 



























l.OE-01 1.0E+00 1.0E+01 I.0E+02 I.0E+03 I.0E-04 

L/E (km/OeV) 



1489.2 days * 25 
1.0E+03 



pu >0.4 GeV 



-*-*■ 
■+ 

+ -*■ 








; \ 
+ 

- + 

■I 

i 


Lu/Ev ++ 4, 


7^ 



l.OE-01 LOE+00 LOE+01 1.0E+02 1.0E+03 l.OE+04 
L/E (km/GeV) 



Figure 25: The ratio of (L ll /E ll ) sc/(L l j,/E ll ) na ii for 14892 
live days (10 SK live days). 



Figure 27: L^/Ey distribution without oscillation for 37230 
live days (25 SK live days). 



2.4-3. L^/E„ } sk distribution for the oscillation 

Instead of E v which is correctly sampled from the 
corresponding probability functions, let us utilize 
E v ,sk which is obtained from the "approximate" 
formula (Eq.(6)) in the preceding paper[l|). 

We express E v described in Eq.(6) utilized by 
Super-Kamiokande Collaboration as E Vt sK to dis- 
criminate from our E v obtained in the stochastic 
manner correctly. 

In Figure [3D1 we give L^/E v ^sk distribution for 
14892 live days (10 SK live days), comparing with 
Lfj,/ 1 E v distribution. It is understood from the com- 
parison that there is no significant difference be- 
tween L^/Ev^sk distribution and L^/E v one in a 
logarithmic scale. This fact tells us that the " aprox- 
imate" formula for E v used by Super-Kamiokande 
Collaboration does not produce so significant dif- 
ference in the logarithmic scale, which may be acci- 
dental. Although the approximated foumula is not 
suitable for the treatment of stochastic quantities 
(see discussion in 3.3 [l|), the result is understand- 
able from Figure 14 in the preceding paper [lj, be- 
cause there is no significant difference between the 
real distribution (correlation) and the " aproximate" 



formula apparently. Also, we can conclude that we 
do not find any hole corresponding to the maximum 
oscillation in L^/E^sk distributions. The reason 
why the Figure can not show such dip structure 
as shown in Figures [3] and [SJ comes from the situa- 
tion that the role of L„ is much more crucial than 
that of E v in the L/E analysis. Namely, L v cannot 
be replaced by L M at all. 

In Figure 1301 we give the comparison of L^/E^ 
distribution with L^/E u ^sk one. 

The apparent small difference between L^j ' E v ^k 
distribution and L^/E v one in Figure [BUI may come 
from that plays an effective role in comparison 
with E v or E^sk, in spite of the situation that 
there are non-negligible differences between E v or 
Eu.sk (see Figure 19 of the preceding paper[l|). 

In Figure 1311 we compare L^jEy distribution 
with Lfj,/ ' E v one. The pretty overlapping between 
them in a logarithmic scale denotes that L M play an 
important role while the energies concerned only 
play the secondary role. The similar situation is 
expected in L v . In Figure [351 we compare L v jE^ 
distribution with L u /E u one. It is clear from the 
figure that we find the first maximum oscillations 
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Figure 28: L^/E^ distribution with oscillation for 1489.2 
live days (one SK live day). 
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Figure 29: L^/Ev distribution with oscillation for 37230 live 
days (25 SK live days). 
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Figure 30: L^j E V ^$K distribution in comparison with 
L^/Ev distribution with oscillation for 14892 day (10 SK 
live days). 



in both distributions on nearly same locations. An- 
other clearer situation is found in the comparison of 
both distributions expressed in a linear scale shown 
as Figures 171 and IT5l It become clear from the com- 
parison of Figures [301 ED an d EH that the flight 
length, either L u or L^, plays a decisively impor- 
tant role in any L/E distribution, compared with 
the energies concerned, either E v or E^ , as it should 
be. 



3. Comparison of L/E Distribution in the 
Super-Kamiokande Experiment with our 
Results 

In our classification, L/E distribution by Super- 
Kamiokande Collaboration should be compared 
with our L^/E^ distribution because they mea- 
sure L M and Efj, directly, or it should be compared 
with our L^/Ev distribution because they get E v 
through the transform from E^. However, they as- 
sert that they measure L v on the SK assumption on 
the direction. Consequently, we compare here their 
L/E distribution with our Ly/E^ one at first. In 
Figure 1531 we compare our L v j ' E v distribution with 
their single ring muon events among Fully Con- 
tained Events should be comparecH as correspond- 
ing ones. 

There are two important matters to be examined 
in the L/E distribution related to the shapes be- 
tween ours and theirs which can be discussed with- 
out entering the details for technical and experi- 
mental condition or criteria around their experi- 
ments. The first one is related to the location and 
its shape for the first maximum oscillation. And the 
second one is related to the location which give the 
maximum frequency of the events concerned. In the 
first one, we can observe the first maximum oscilla- 
tion at L v /E v = 515 km/GeV (see Figures H to UHl 
exactly in our computer numerical experiment un- 
der the oscillation parameters obtained by Super- 
Kamiokande Collaboration. Furthermore, we can 
observe clearly the second, the third or more higher 
order maximum oscillations at the anticipated lo- 
cations (see Eq.(3) and Figure [7]) in our numerical 
experiment. Also, the shapes of those maximum 
oscillations are rather sharp which comes from the 



4 We read out Fully Contained Events among total events 
from Super-Kamiokande Collaboration QQ- Since we are 
only interested in single ring muon events, events with the 
highest quality, excluding Partially Contained Events for our 
analysis. 
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Figure 31: Comparison between L^/Eu distribution and 
L/J./E/J, distribution with oscillation for 37230 days (25 SK 
live days). 
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Figure 32: Comparison between L y /E v distribution and 
Liz/E^ distribution with oscillation for 37230 days (25 SK 
live days). 
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Figure 33: The comparison of L/E distribution for single- 
ring muon events due to QEL among Fully Contained Events 
with the corresponding one by the Super-Kamiokande Ex- 
periment. 



specified oscillation parameters obtained by Super- 
Kamiokande Collaboration. On the other hand, 
they obtain a broader region for the absence of the 
neutrino events as the result of the first maximum 
oscillation such as 100 < L/E < 800 (km/GeV). 
Such a broader region may contradict the concept 
of the survival probability for a given flavor un- 
der the specification of their oscillation parameters, 
taking account of the results from the analysis of 
the single ring muon events among Fully Contained 
Events, the highest quality events among all events 
to be analyzed by them. 

Now, we examine the remarkable difference be- 
tween ours and theirs as for the locations of the 
maximum frequencies for the events. As shown 
in Figure [33J we g ive it as 1.0 < L u /E u < 1.26 
(km/GeV), while they give 20 < L v /E v < 25 
(km/GeV) which is larger than ours by one order 
of the magnitude. 

Here, at first, in our computer numerical exper- 
iment, we discuss the correlation between L and 
E at the location for the maximum frequency for 
the events and the similar correlation at the cor- 
responding location where Super-Kamiokande Col- 
laboration give their maximum frequency. Next, we 
clarify what happen at the location for the maxi- 
mum frequency for the events in their experiment. 

In the following discussion, we designate neutrino 
events with 135° < 6 V < 180° as vertical like events, 
events with 90° < B v < 135° as horizontal like 
events and events with 88° < Q v < 92° as exclu- 
sively horizontal events, respectively. 

In Figures [3H to [351 we give the correlations 
for the maximum frequency for the events in our 
L v jE v distribution. Also, in Figures 1371 to 1391 we 
give the similar correlations for the location in our 
L v / E v distribution which correspond to the max- 
imum frequency of the events obtained by Super- 
Kamiokande Collaboration. In Figure [331 we gi ve 
the correlation between L v and E v for the interval 
1.0 < L v jE v < 1.26 (km/GeV) at the maximum 
frequency for the events in our computer numeri- 
cal experiment. It is clear from the figure that all 
incident neutrino events have the values of L v less 
than 10 km, corresponding to the vertical-like, hori- 
zontal like and the exclusively horizontal downward 
neutrinos, taking account of the transform between 
Lu{Lfj,) and cos9 v {cos9^), as they must be. All 
these neutrinos cover all incident directions as the 
downward and they correspond to the maximum 
frequency for the events. It is clear from Figure [371 
that the incident neutrino events are concentrated 
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Figure 34: Correlation diagram between L v and E v for 
1.0 < L v /E v < 1.26 (km/GeV) which corresponds to the 
maximum frequency of the neutrino events for L v /E v dis- 
tribution in our computer numerical experiment for 14892 
live days (10 SK live days). 
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Figure 37: Correlation diagram between L v and E v for 20 < 
Lv/Ev < 25 (km/GeV) which corresponds to the maximum 
frequency of the neutrino events for L^/Eu distribution in 
SK experiment for 14892 live days (10 SK live days). 
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Figure 35: Correlation diagram between L v and L M for 
1.0 < Lu/Ev < 1.26 (km/GeV) under the neutrino oscil- 
lation parmeters obtained by Super-Kamiokande Collabora- 
tion for 14892 live days (10 SK live days). 
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Figure 38: Correlation diagram between L v and L^ for 
20 < L v /E v < 25 (km/GcV) under the neutrino oscillation 
parmeters obtained by Super-Kamiokande Collaboration for 
14892 live days (10 SK live days). 
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Figure 36: Correlation diagram between L M and E M for 
1.0 < L v /E v < 1.26 (km/GeV) which corresponds to the 
maximum frequency of the neutrino events for L v /E v dis- 
tribution in our computer numerical experiment for 14892 
live days (10 SK live days). 



Oscillation Am 2 = 2.4x 10 J eV 2 20 < Lv/Ev < 25. 
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Figure 39: Correlation diagram between L M and E M for 20 < 
L v /E v < 25 (km/GcV) which corresponds to the maximum 
frequency of the neutrino events for L v /E v distribution in 
SK experiment for 14892 live days (10 SK live days). 
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Figure 40: Correlation diagram between L M and E v for 
15.8 < L^/Ey < 31.6 (km/GeV) which correspond to the 
maximum frequency of the neutrino events for L^/E v dis- 
tribution in SK experiment for 14892 live days (10 SK live 
days). 
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Figure 41: Correlation diagram between L u and L M for 
15.8 < L^/E v < 31.6 (km/GeV) under the neutrino os- 
cillation parmeters obtained by Super-Kamiokande Collabo- 
ration for 14892 live days (10 SK live days). 
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Figure 42: Correlation diagram between L v and E^ for 
15.8 < L^/Ey < 31.6 (km/GeV) which correspond to the 
maximum frequency of the neutrino events for L^/E v dis- 
tribution in SK experiment for 14892 live days (10 SK live 
days). 



into the interval of from 20 to 100 km for the value 
of L u which correspond to the horizontal like events 
and the exclusively horizontal events, not including 
vertical like events, taking the account of the rela- 
tion between L V (L^ and cos6 v {cos0^). It is quite 
natural that events number there is smaller than 
that for the maximum frequency which includes the 
vertical like events. 

In Figure 1351 we give the correlation between L v 
and for the same intervals as in Figure [34] It 
is clear from the figure that the majority of the 
events is concentrated into the squared region with 
L u < 10 km and < 10 km. This denotes that 
the downward incident neutrinos produce muons in 
the forward direction irrespective of scattering an- 
gles. At the same time, it should be noticed from 
the figure that the non-negligible parts of the down- 
ward incident neutrino events produce the upward 
muons (^1000 to 10000 km) due to the backscat- 
tering as well as the azimuthal angle effect due to 
QEL for both horizontal like and exclusively hori- 
zontal events which is clearly shown in Figure 1361 
too. These upward muons may be surely identi- 
fied as the products of the upward neutrinos in the 
analysis performed by the Super-Kamiokande Col- 
laboration. 

Furthermore, from the comparisons of Figure [35] 
with Figure [33J and of Figure [35] with Figure 1391 
it is easily understood that exclusively horizontal 
neutrinos which occupy the majority in Figures 1351 
to [39] are more influenced by the effects of both 
the backscattering and the azimuthal angle in QEL, 
compared with the cases in Figures 1331 and 1361 This 
denote exclusively horizontal like neutrino (down- 
ward) produce upward muons (backward direction) 
with higher probability compared with both ver- 
tical like events and horizontal like events. Thus, 
from the comparison of Figures [3J] to [33] with Fig- 
ures [37] to [SHI it is concluded that there is no contra- 
diction for the interpretation of all figures between 
the maximum frequency for the event in our L v /E v 
distribution and the corresponding distribution at 
the location where Super-Kamiokande Collabora- 
tion give their maximum frequebcy. 

Now, we examine the reliability of the maxi- 
mum frequency of the events obtained by Super- 
Kamiokande collaboration as shown in Figure 1331 
They assert that they measure the directions of 
the incident neutrinos by measuring those of muons 
under the SK assumption on the direction. How- 
ever, what they measure really are the directions 
of the muons, but not those of the corresponding 
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neutrinos due to their neutrality. Consequently, 
here, we examine the L^/E^ distribution in de- 
tail, not L v /E v distribution for checking the exper- 
imental data obtained by Super-Kamiokande Col- 
laboration. Strictly speaking, they measure 
and E^, not and E v . However, as they trans- 
form the original E^ to E v (Eq.(7) in Parti |lj), 
we interpret they "measure" L M and E v . Thus, 
we examine our L^/ E v distribution for the inter- 
val 15.8 < L v /E v < 31.6 (km/GeV) where Super- 
Kamiokande Collaboration give their maximum fre- 
quency of the events. 

In Figures HO] to EH we give our correlations for 
the events at the location 15.8 < L^/E^ < 31.6 
(km/GeV) where Super-Kamiokande Collaboration 
give their maximum frequency for the events. In 
Figure 001 we give the correlation between and 
E u . In Figure SU we give the correlation between 
L v and L M . In Figure |42j we give the correlation 
between L v and E^. 

It is understood from these figures that here, 
neutrino events produce exclusively the down- 
ward muons which consist of horizontal like events 
and exclusively horizontal events, but not vertical 
like events, taking account of the transform from 
L U (L^) to cosOy (cosOfi). Also, it is easily seen that 
such the muons are produce by the parent neutri- 
nos whose zenith angle distribute over downward to 
upward widely. This fact shows without doubt that 
one cannot decide the direction of the incident neu- 
trino even the maximum frequency for the events 
from the measurement of the produced muons. 

Here, we comment to the recent work on 
L/E analysis by Dufours^Q, a member of Super- 
Kamiokande Collaboration^]. In her paper, she has 
carried out the Monte Carlo simulation with neu- 
trino oscillation around L/E analysis and has ob- 
tained a beautiful agreement between the exper- 
imental data and her Monte Carlo results. This 
seems to be the first Monte Carlo simulation 
with oscillation in Super-Kamiokande Collabora- 
tion, since before this work, Super-Kamiokande 
Collaboration have been comparing their exper- 
imental results with their Monte Carlo simula- 
tion without oscillation and have estimated neu- 



5 One may take a notice that it is not appropriate to cite 
Ph.D due to their nature of "unpublished". However, many 
Ph.D thesis around Super-Kamiokande have been publised 
and there are no any contradictions between their Ph.D the- 
ses and SK papers and the detailed descriptions are exclu- 
sively found in these Ph.D theses. 
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Figure 43: The excluded region for L/E analysis of FC 
single-ring events by Super-Kamiokande. 



trino oscillation parameters from the difference be- 
tween them. In order to keep the consistency with 
the usual Monte Carlo simulations without oscil- 
lation performed by Super-Kamiokande Collabora- 
tion, she must have been carried out her Monte 
Carlo simulation with oscillation under the SK as- 
sumption on the direction, because the SK assump- 
tion on the direction is the cornerstone throughout 
their analysis on neutrino oscillation. It is too clear 
that the results obtained by us contradict her result, 
even if considering the difference that we examine 
the Fully Contained Events only, while she has ex- 
amined the Partially Contained Events in addition 
to the Fully Contained Events. Furthermore, it may 
be unnatural that she has obtained extremely beau- 
tiful agreement with the experimental data. 

Finally, we examine the data selection procedure 
made by Super-Kamiokande (hereafter called as SK 
data selection procedure) which is imposed upon 
their Fully Contained Events in the single ring muon 
events. They introduce such a procedure to ex- 
clude ambiguity mainly coming from horizontal like 
events. They exclude single ring muon events as 
Fully Contained Events which exist within the re- 
gion described in Figure 2(a) in their paper We 
reproduce it in Figure l43l The region enclosed by 
two lines is the excluded region in their analysis. 
Here, it should be emphasized that we need not ex- 
clude the horizontal like neutrino events at all in 
our L/E analysis against the SK analysis, because 
we have not any experimental errors even in our 
horizontal-like neutrino events due to the nature of 
our computer numerical experiment. However, it is 
not vain to examine whether the SK data selection 
procedure is appropriate or not for our analysis in 
our computer numerical experiment. In Figure l44l 
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Figure 44: The relation between the excluded region and 
the correlation between L v and E v for 10 SK live days. 
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Figure 45: Lv/E u distributions with and without the re- 
striction imposed by Super-Kamiokande for 10 SK live days. 
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Figure 46: The relation between the excluded region and 
the correlation between L M and for 10 SK live days. 
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Figure 47: L^/E^ distributions with and without the re- 
striction imposed by Super-Kamiokande for 10 SK live days. 



we give correlation diagram between L v and E v 
with oscillation together with the lines for exclu- 
sion given by Supper-Kamiokande Collaboration. 
The lines for exclusion in Figure |44] are transformed 
from the lines in Figure S3] The region enclosed 
by two lines in Figure [U are the excluded region. 
Namely, the events concerned within the excluded 
region should be subtracted from the original distri- 
bution. Thus, L v jE v distribution of the resultant 
events after subtraction procedure (with restriction 
in the figure) is shown in Figure 25] together that 
without subtraction (no restriction in the figure). 
It is clear from the figure that L v jEv distribution 
with the SK data selection procedure keeps still the 
characteristics of a series of the maximum oscilla- 
tions in spite of the decrease of the events in smaller 
value of L v / E v . Here, we show L v j ' E v distribution 
under the SK data selection procedure in a linear 
scale to make their physical image clearer. It is 
needless to say that we need not introduce the SK 
data selection procedure into our computer numer- 
ical experiment due to the "complete experiment" 
by the definition, but even if we introduce it, Fig- 
ureSS]shows that the essential characteristics of our 
L v / E v distribution is never changed. 

Similarly we examine the case of L^/E^ distri- 
bution. Figure 05] for L^/E^ distribution corre- 
sponds to Figure |4"4"1 for L u /E u distribution. Also in 
Figure U7] L^/E^ distribution corresponds to Fig- 
ure 05] for L V IE V distribution. 

It is clear from Figure [47] that the L^/E^ distri- 
bution with the SK data selection procedure do not 
show anything like the maximum oscillation in the 
same way as the original L^/E^ distribution. 

Summarized from Figures 03] to 07] we could not 
extract the neutrino oscillation parameters from 
L^/E^ distribution, even if we apply the SK data 
selection procedure to the original L^/E^ distribu- 
tion, while we can keep the essential character of 
the maximum oscillations in L v / E u analysis, even 
if we apply to the SK data selection procedure to 
the original Lv/E v distribution. Finally, we should 
say again that we need not introduce the SK data 
selection procedure into our computer numerical ex- 
periment due to the nature of no error experiment. 

4. Conclusion 

The determination of the neutrino oscillation pa- 
rameters is entirely based on the survival proba- 
bility for a given flavor. Consequently, it is in- 
evitable to decide L„ and E v as precisely as pos- 
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sible, when one wants specify neutrino oscillation 
parameters. However, in cosmic ray experiments, 
one may not measure both L v and E v due to their 
neutralities, in addition to because of the incapa- 
bility for the determination of the directions of the 
incident neutrinos due to the essential nature of cos- 
mic ray beams. Therefore, cosmic ray physicists are 
forced to assume the direction of incident neutrinos 
a priori, when they do not carry out computer nu- 
merical experiment as their second experiment. In 
the case of Super-Kamiokande Collaboration, they 
assume that the direction of incident neutrino is 
approximately the same as that of emitted lepton 
(the SK assumption on the direction). In the pre- 
ceding paper we have verified that the SK as- 
sumption on the direction does not hold even if 
approximately. The essential conclusion obtained 
by the present paper is that in principle one may 
not specify the neutrino oscillation parameters from 
the cosmic ray experiments due to the unknown 
directions of the incident neutrinos. Our verifica- 
tion that Super-Kamiokande Collaboration cannot 
specify the neutrino oscillation parameters through 
their L/E analysis at least, and consequently, our 
conclusion can be applied to any type of experiment 
of cosmic ray physics where the direction of the neu- 
trino, in principle, cannot be determined. Our con- 
clusion tells us that only accelerator physics can 
specify the neutrino oscillation parameters reliably, 
if the neutrino oscillation really exists. 
Deduction of our conclusion is as follows: 

(1) There are much uncertainty factors in cos- 
mic ray physics, compared with the acceler- 
ator physics due to the original nature of cos- 
mic ray. Consequently, in spite of such a diffi- 
culty, if one wants to carry out the experiment 
with high precision on neutrino oscillation, we 
should focus the simplest and clearest "target" 
by which one get high quality information on 
neutrino oscillation. With such a motivation, 
we choose the single ring muon events due to 
QEL which they occur inside the detector and 
terminate inside the detector (Fully Contained 
Events). Here, the kind of neutrino concerned 
is clear (electron neutrino or muon neutrino). 
The energy of emitted lepton and its direction 
can be estimated reliably ( from the standpoint 
of Super-Kamiokande at least). The circum- 
stance around our computer numerical experi- 
ment is modeled after real Super-Kamiokande 
experiment in essential points. We have ana- 



lyze the single ring muon events due to QEL as 
Fully Contained Events obtained by our com- 
puter numerical experiment. 

(2) We have carefully and in detail examined the 
validity of the SK assumption on the direc- 
tion which is the " cornerstone stone" for their 
analysis around neutrino oscillation. As the 
result, we have clarified that the SK assump- 
tion on the direction does not hold even if ap- 
proximately. Also, we examine the validity of 
the unique relation between E v and E^ ex- 
pressed in the polynomial obtained by Super- 
Kamiokande Collaboration and we have clar- 
ified that the unique relation between them 
does not hold. These two improper treatments 
originate from the situation that they do not 
consider characteristics of the stochastic pro- 
cesses concerned seriously. However, the unre- 
liability on the directions of the incident neu- 
trinos influences final result in a fatal manner, 
while unreliability on the energy estimation 
does not provide the significant errors com- 
pared with the former. The concrete sum- 
maries are given in (3). 

(3) Due to the nature of the computer numerical 
experiment, assuming neutrino oscillation pa- 
rameters obtained by Super-Kamiokande Col- 
laboration, we carry out all possible combina- 
tion of L/E analysis, namely, L v /E v analysis, 
LyjE^ analysis, L^/E v analysis and L^/E^ 
analysis, based on the survival probability for 
a given flavor whose variables are L„and E v . 
Among four L/E analyses, only L v /E v analy- 
sis has reproduced the existence of the maxi- 
mum oscillations, not only the first maximum 
oscillation but also the second, the third, the 
fourth and so on. The confirmation of a se- 
ries of the maximum oscillations, such as ,the 
first, the second, the third and so on in L v /E v 
analysis shows that our computer numerical 
experiments have been carried out in a correct 
manner. The L v /E^ analysis has reproduced 
the first maximum oscillation roughly, but can- 
not reproduce the maximum oscillation after 
the second. Both L^/Ey and L^/E^ analy- 
ses cannot reproduce any characteristics of the 
maximum oscillation at all. Notice that Super- 
Kamiokande Collaboration have carried out ei- 
ther L^/Ey analysis or L^/E^ analysis, neither 
L v /E v analysis nor L v / E^ analysis. Combined 
with the item (1), these facts tell us that the 
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decisive variable in the survival probability for 
neutrino oscillation is L u but neither L^, nor 
E v , nor Ep. Thus, our verification that the 
SK assumption on the direction does not hold 
even approximately requests urgently Super- 
Kamiokande Collaboration to re-analyse their 
results around the zenith angle distributions 
for neutrino events which have been regarded 
as the establishment of the existence of the 
neutrino oscillation, because their analysis on 
neutrino oscillation (atmospheric neutrino) is 
entirely based on the survival probability and 
the SK assumption on the direction. 

Finally, we would like to emphasize the im- 
portance of the cosmic ray study in order to 
avoid any misunderstandings. This character- 
istics of the cosmic ray study never make it lose 
its raison-d'etre. The main role of cosmic ray 
physics is to grasp qualitatively the essential 
of something like new. Up to now, cosmic ray 
study has been contributing to find new indi- 
cations in fundamental physics and from now 
on, it will be continue to fulfill its role. 
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